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Abstract 

A supersymmetric method for the construction of so-called condition- 
ally exactly solvable quantum systems is reviewed and extended to classical 
stochastic dynamical systems characterized by a Fokker-Planck equation with 
drift. A class of drift-potentials on the real line as well as on the half line 
is constructed for which the associated Fokker-Planck equation can be solved 
exactly. Explicit drift potentials, which describe mono-, bi-, meta-or unstable 
systems, are constructed and their decay rates and modes are given in closed 
form. 

1 Introduction 

Many physical phenomena of nature are characterized by some basic differential 
equations. For example, quantum-mechanical phenomena are described by Schro- 
dinger's equation, which dictates the dynamics of some quantum system represented 
by a Hamilton operator. One is therefore primarily interested in finding all eigenval- 
ues and eigenstates of such Hamiltonians. As a consequence, finding a large class of, 
in this sense, analytically exactly solvable quantum systems is an important goal and 
this search has already been initiated by Schrodinger using the so-called factoriza- 
tion method [|1], |2[. This factorization method is basically equivalent to Darboux's 
method already invented 1882 and applied to Sturm-Liouville problems. Re- 
cently, these methods have attracted new attention [|J, in particular in connection 
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with supersymmetric (SUSY) quantum mechanics ||. For a recent approach and 
further references see, for example, [|, [7], || . 

Somehow unnoticed by those working on the quantum mechanical problem the 
same basic ideas have been utilized in constructing exactly solvable classical stochas- 



tic systems described by a Fokker-Planck equation || [TI| [TT|] - In fact, the Fokker- 
Planck equation can be put into the form of an imaginary-time Schrodinger equation 
exhibiting SUSY. Therefore, exactly solvable SUSY quantum Hamiltonians also pro- 
vide exactly solvable Fokker-Planck-type systems. Whereas until now these methods 
have been applied to mono- and bistable systems || [LO], [TTH , they correspond to a 
situation with unbroken SUSY, we present here also examples with broken SUSY 
giving rise to meta- or unstable drift potentials. 

In the next two sections we briefly review the basic concepts of SUSY quan- 
tum mechanics and its close connection with the Fokker-Planck equation ||, |P2"f . 
Section 4 reviews our recent approach || [7], § towards a construction of the most 
general class of SUSY partners of a given quantum mechanical Hamiltonian. This 
approach is applied to the Fokker-Planck equation and allows to construct new drift 
potentials with known decay rates and modes. Two examples are discussed in some 
detail. The first one is related to the linear harmonic oscillator having unbroken 
SUSY and leads to bistable double-well or monostable single-well drift potentials 
already found by Hongler and Zheng || . The second example is the class related to 
the radial harmonic oscillator with broken SUSY. Here the drift potential is either 
metastable or unstable. 



2 Supersymmetric quantum mechanics 

To begin with, let us briefly review the basics of SUSY quantum mechanics or, to be 
more precise, Witten's model of SUSY quantum mechanics |J. This model consists 
of a pair of standard Schrodinger Hamiltonians (units are such that h = m = 1) 

H± = -~ + V ± (x) , V ± (x)= 1 -(w 2 (x)±W'(x)) . (1) 

Evidently, this pair is completely characterized by the so-called SUSY potential W, 
which is assumed to be real- valued and an at least once differentiable function inside 
the configuration space. This configuration space will be either the real line R or the 
positive half line M + . Thus the above Hamiltonians act on a suitable linear space 
of square-integrable functions over R and 1R + , respectively. With the help of the 
supercharge operators 

the above pair of SUSY-partner Hamiltonians factorizes, H + = AA^ > 0, H = 
A^A > 0, and obviously obeys the intertwining relations AH_ = H + A and H_A^ = 
A J[ H + . As a consequence H + and H_ are essentially isospectral. In other words, 
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their strictly positive energy eigenvalues coincide and the corresponding eigenstates 
are related via the supercharges. However, there may exist an additional vanishing 
eigenvalue for one of these Hamiltonians. In this case SUSY is said to be unbroken 
and by standard convention |5j this additional eigenvalue is assumed to belong to 
H_. If such a state does not exist SUSY is said to be broken. To summarize, for 
unbroken SUSY we have 

E =0, E~ +1 = E+>0, ^(x)=i/>o(0)exp{-J*dzW(z)j, (3) 
= {EtY l ' 2 Am , V+ = (E- +1 )-^A^- +1 , 
whereas for broken SUSY these relations read 

E~=E+>0, C = (£n + r 1/2 ^, 1# = {E-y^A^ . (4) 

Here we have denoted the eigenfunctions and eigenvalues of H± by tp^ and E^, 
respectively. That is, 

H±if>£ = E±H , 71 = 0,1,2,... , (5) 

and we also note that throughout this paper we will consider, without loss of gen- 
erality, quantum systems with a purely discrete spectrum. 

We conclude this section by noting that for broken as well as unbroken SUSY 
one can obtain the complete spectral information of one Hamiltonian, say H_, if the 
eigenvalues and eigenstates of the corresponding partner, here H + , and the SUSY 
potential W are known. 



3 Classical stochastic dynamics 

As mentioned in the Introduction we are also interested in classical systems with a 
stochastic dynamics governed by the Fokker-Planck equation 



— m t (x,x ) 



1 d 2 



d ( 

> r7r 2' mi ( x > x °) + -^(U'{x)m t (x,x 



(6) 



Here m t (x, x ) denotes the transition-probability density of a macroscopic degree 
of freedom to be found at time t at position x if it initially has been at Xq, i.e. 
mo(x, xq) = 6(x — x ). This degree of freedom is subjected to an external force 
characterized by the real-valued drift potential U and to a stochastic random force 
(white noise) resulting in the diffusive term on the right-hand side of (^|) with diffu- 
sion constant set equal to 1/2. Making the ansatz |12[ 



m t (x,x ) = e u{x) K t (x,x ) 



leads to 



1 d 2 



2 dx 2 



12/ 



U"{x) 



K t (x,x ) 



(7) 



(8) 
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which may be interpreted as an imaginary-time Schrodinger equation for the SUSY 
Hamiltonian if_ with a SUSY potential given by the first derivative of the drift 
potential, W = U'. Hence, the desired transition-probability density is given via the 
Euclidean propagator for H_: 

m t (x,x ) = exp{U(x ) - U(x)}(x\ exp{-tH^}\x ) . (9) 

It is obvious that the decay modes and decay rates of this classical stochastic dy- 
namical system are related to the eigenfunctions and eigenvalues of To be more 
explicit, let us consider the cases of unbroken and broken SUSY separately. 

For unbroken SUSY the ground-state energy of H_ vanishes and as a consequence 
there exists a stationary, i.e. time-independent, probability distribution. Noting that 
ipo(x) = iPq{xq) exp{U(x ) — U(x)} the transition-probability density can be put 
into the form 

oo 

m t (x,x ) = {x)] 2 + exv{U{xv) - U{x)}^e- tE " i)-{x)i)-{x Q ) , (10) 

71=1 

which clearly shows that the strictly positive eigenvalues of H_ and the associated 
eigenfunctions characterize the decay rates and decay modes of the system. In 
addition the ground-state wavefunction represents the stationary distribution. That 
is, a stable system is characterized by an unbroken SUSY. 

In the case of broken SUSY, by definition, if_ does have strictly positive eigen- 
values only and cannot lead to a stationary distribution. In other words, broken 
SUSY is related to unstable systems. Here the transition probability density has 
the form 

oo 

m t (x,x ) = exp{U{x ) - U{x)}^2e' tE "^j-{x)ij-{xo) . (11) 

71=0 

Finally, we mention that in both cases we may invert the drift potential, U — > 
—U, which amounts in replacing H_ by H + . Hence, due to the presence of SUSY 
the inverted drift potential has the same decay rates as the original one. Only in 
the case of unbroken SUSY a stable system becomes unstable upon inversion. This 
is expressed in the fact that the vanishing eigenvalue Eq = is missing in the spec- 
trum of H + . The decay modes of the original and the inverted drift potential are 
also clearly related by the SUSY transformations @ and respectively. 

4 Designing exactly solvable models 

From the discussion in the previous two sections it is clear that once the spectral 
data, that is, the eigenvalues and eigenfunctions, of H + are known we also know 
these data for the corresponding superpartner if_ and in turn the decay rates and 
modes for the Fokker-Planck equation with drift potential U(x) = dzW(z). In 
order to construct new exactly solvable systems we will search for the most general 
class of superpartners associated with a given exactly solvable quantum Hamiltonian 
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H + . Such a construction method has recently been given [7|, |Sj and we briefly 
review it here. 

In order to find the most general class of SUSY partners of a given SUSY Hamil- 
tonian we make the following ansatz 

W(x) = *(x) + , (12) 

where $ is a known SUSY potential belonging to some shape-invariant (i.e. exactly 
solvable) potential. See, for example, Table 5.1. in ref. f|. If we now assume that u 
is a solution of 

u"(x) + 2<$>(x)u'(x) -bu(x) =0 , beR, (13) 

we find that 

V + (x) = U\x) + U'(x) + b -. (14) 

By construction V + is up to the additive constant b/2 shape-invariant and, therefore, 
the eigenvalues and eigenf unctions of H + are exactly known. The corresponding 
partner potential can be put into the form 

Tr / x 1 *.n \ 1 \ u'(x) / x , . u'(x)\ b 

= 2 * W ~ 2 * W + iJ) V * W + W)) - 2 ' <15) 

which, in general, is not shape-invariant and thus a new exactly solvable quantum 
mechanical potential. Similarly we may construct also a new drift potential 

U(x)= I &z<$>(z) + \ogu(x) (16) 

Jx 

for which the associated decay rates and modes are known exactly. 

Let us note here that we cannot take an arbitrary solution of ([13|) for the con- 
struction of new quantum-mechanical and drift potentials. In order to circumvent 
domain questions of the operators involved we will take only those solutions of ( |T3"D 
into account which are strictly positive. Hence, no singularities inside the configura- 
tion space will appear in V- and U. For this reason the new Schrodinger potentials 
U_ obtained in this way have been called conditionally exactly solvable ||. In the 
following two subsections we will present two examples with unbroken and broken 
SUSY, respectively. For further details and examples see ||. 

4.1 A system with unbroken SUSY 

The first example we are going to present is characterized by a linear SUSY potential, 
$(x) = x, with configuration space given by the real line, x G R. This SUSY 
potential gives rise to the harmonic-oscillator potential 

V+(x) = \(x 2 + b + l) (17) 
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and the corresponding eigenvalues and eigenfunctions of H + are 

E+ = n + b/2 + 1 , V« 0*0 = [v^2 n n!] ~ 1/2 # n (z) exp{-x 2 /2} . (18) 

Here H n denotes a Hermite polynomial of order n. In this case the general solution 
of ( pT3|) can be expressed in terms of confluent hypergeometric functions, 

u{x) = X F X (-f , |, -x 2 ) + /3x xFx (2f*, |, -x 2 ) (19) 

and is strictly positive for b > -2 and < 2r(| + l)/r(^), cf. ref. 0, g. The 
corresponding partner potential is given by 



VJx) = -x 2 h —H- 

K ' 2 2 



u'(x) 
2x + —Hr 



111 X I 



(20) 



and plots of it for various values of the parameters h and (3 can be found in Figure 1 
of [^] and Figures 1 and 2 of ||. Here we only note that for large \x\ this potential 
becomes asymptotically that of a harmonic oscillator. For values of b close to the 
lower bound —2 this potential exhibits two double-wells near the origin. Whereas 
for larger values of b these double wells merge to a single well at the origin. For 
(3 = the potential V_ is symmetric about x = 0. This symmetry is broken for non- 
vanishing (3. Noting that SUSY is unbroken for all allowed values of the parameters, 
the spectral properties of the corresponding Hamiltonian H_ are easily obtained 
from (0): 

E~ = , E- +1 = Ei = n + b/2 + 1 , ^~(x) = exp{-x 2 /2} , 

VWi(z) = eXP{ ' XV2} ^ + H n (x) ^ 

[ v ^F2«+ 1 n!(n + 6/2 + l)] 7 V 

(21) 

These results also allow a complete study of the Fokker-Planck equation for the 
drift potential 

U{x) = l -x 2 + logw(x) . (22) 

In fact, this drift potential is stable and has a stationary distribution given by 
[iJjq(x)} 2 . The decay rates and decay modes are given explicitly by E~ +1 and VvT+i m 
(pip. Plots of this drift potential for various values of the parameters can be found 
in Figures 1 of . Here we again briefly mention that for b close to its lower limit U 
has the shape of a bistable (double-well) potential being symmetric about the origin 
for (3 = 0. For larger values of b this drift potential develops also a stable single- well 
shape. 

To conclude this subsection let us mention that the results presented here are 
not new. They have first been derived by Hongler and Zheng in 1982, which was 
even before the work of Mielnik [Q who was searching for more general factoriza- 
tions of the harmonic oscillator Hamiltonian. His results can be viewed as special 
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cases of those presented here. The advantage of the present approach is that it can 
be applied not only to the harmonic oscillator-like systems but to all exactly solv- 
able ones || . In particular, there exist also shape-invariant quantum systems which 
allow for a broken SUSY ||. That is, one may be able to design drift potentials 
which are metastable and can be solved exactly. To our knowledge, such exactly 
solvable systems are not available so far. Except, of course, the rather trivial case of 
a piecewise linear drift potential. In the next subsection we are going to present an 
exactly solvable metastable drift potential, which is related to the radial harmonic 
oscillator Hamiltonian. 



4.2 A system with broken SUSY 

In this subsection we will consider the case where the SUSY potential is given by 

= x + - , 7 > , (23) 

x 

where the condition put on the parameter 7 leads to a SUSY potential $ which 
characterizes a shape-invariant SUSY pair of radial harmonic oscillator-like Hamil- 
tonians with broken SUSY. Due to the condition flT3] ) even the SUSY potential ( |T2"D 
with above $ gives rise to the radial harmonic oscillator potential 

which is exactly solvable and leads to the following spectral properties of H + : 

2n\ 



£+ = 2n + 2 7 + l + ^, ^(x) 



r(n + 7 + l/2) 



1/2 

x 7 e -* 2 /2 Lty-^ix 2 ) . (25) 



Here Ln denotes an associated Laguerre polynomial of degree n and index v. As in 
the previous case the general solution of (|13|) can be expressed in terms of confluent 
hypergeometric functions. Here, however, we are only interested in cases with broken 
SUSY and with this constraint the most general solution reads 

u(x) = iFxC-J, 7 + |, -x 2 ) = e- 2 1 F 1 ( 7 + *±*, 7 + f , x 2 ) . (26) 

This solution will be strictly positive if b > —47 — 2. The corresponding SUSY 
partner potential reads 

V - [X) 2 + 2x 2 +7 2 + u(x) V x u(x) 



and is plotted, for example, as Figure 5 in ref. ||. The eigenvalues of the Schrodinger 
Hamiltonian if_ for this potential are identically to those of H + given in (|2"5j). The 
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eigenfunctions can also be obtained from those in ( P5|) via the SUSY transformation 
(|) and read 



2n\ 



{n + 1 + \ + \)T(n + 1 + l/2) 



1/2 

x^ +1 e- x2 ' 2 (l^ +1 ' 2 \x 2 ) + U ' {X) 



2xu[x) 
(28) 

Again, as in the previous case we can construct a drift potential with decay rates 
and modes given by the eigenvalues and eigenfunctions of This drift potential 
is explicitly given by 

U(x) = ix 2 +7logx+logu(x) = -|a; 2 +7logx+log[iFi(7 + ^,7 + i,x 2 )] (29) 

and because of broken SUSY does not have a stationary distribution. To be more ex- 
plicit, for small x > it has a logarithmic "hole" at the origin, i.e. U(x) ~ — 7I logx| 
for 1 < 1, whereas for large x — > 00 it becomes asymptotically that of a harmonic 
oscillator, i.e. U(x) ~ x 2 /2 for x ^> 1. In addition to that, for values of b close to its 
lower limit —47 — 2 this drift potential exhibits a local minimum. In other words, 
it is metastable. In fact, the smallest decay rate is given by Eq =27 + l + 6/2\0 
for b \ —47 — 2. For larger values of b this local minimum disappears and the drift 
potential ( |2"9"| ) becomes unstable. A typical plot of this potential is given in Figure 
1 where we have shown ( f2~9"| ) for 7 = 1 and b e (—6, —5). Note that in Figure 1 we 
have plotted U(x) versus exp{— x}. 



5 Final remarks 

In this paper we have extended a recent approach to the construction of drift po- 
tentials for which the associated Fokker-Planck equation can be solved exactly. As 
starting point we have chosen a SUSY potential $ which generates a pair of shape- 
invariant quantum mechanical potentials leading to exactly solvable quantum Hamil- 
tonians. This SUSY potential can be perturbed (by adding u'/u) in such a way that 
one of these Hamiltonians, H + , remains in the class of shape- invariant exactly solv- 
able Hamilton operators. The partner Hamiltonian if_, however, is a new one and 
due to SUSY its spectral properties can be obtained from H + in a straightforward 
way. Using the close relations between SUSY quantum mechanics and the Fokker- 
Planck equation we have used these results to find also new drift potentials with 
exactly known decay rates and modes. In this paper we have constructed stable as 
well as unstable drift potentials associated with unbroken and broken SUSY, respec- 
tively. The examples for unbroken SUSY are actually not new and have already been 
studied by Hongler and Zheng ||. However, the present results for broken SUSY 
leading, in particular, to metastable drift potentials are new. To our knowledge, 
these are the first (on the positive half line) analytical metastable drift potentials 
for which the decay rates and modes can be given in closed analytical form. Besides 
some practical applications these metastable potentials can also serve as a testing 
ground for approximation methods. Recall that for metastable drift potentials the 



8 



Figure 1: A family of drift potentials (|29|) for 7 = 1 showing a transition from 
metastability to instability with increasing b. The corresponding decay rates and 
decay modes are known in closed form. Note that we have plotted U(x) versus 
exp{— x}. 



fluctuation operator for the classical "bounce" solution has a negative eigenvalue 
and, therefore, leads to serious singularities within a saddle-point approximation 
T3fl . The standard treatment in such unstable situations goes back to Langer [14 



and is based on some analytical continuation techniques. These formal treatments 
can now be applied to and (for the first time) tested with the examples presented 
in Section 4.2. 
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